Claims: 

1. A method for registering images, comprising the steps of: 
receiving a first image and a second image, said first image and said second 

image including three-dimensional data sets; 

globally transforming said second image via a rigid transformation techniqu; 
locally transforming said first image via an iterative motion tracking 
technique; and ^ 

outputting a registered iniage. 

\ 

2. The method of claim 1, wherein said globally transforming includes the steps 
of: 

constructing a first pyramid for ^aid first image; 

constructing a second pyramid foresaid second image; 
15 determining a set of optimal scale fkctors; 

performing a uniform scaling on saicksecond image in accordance with said set 
of optimal scale factors to obtain a rescaled second image; 

determining a set of optimal translation and rotation parameters; and 

performing a transformation on said rescai^d second image based on said set 
20 of optimal translation and rotation. 

3. The method of claim 2, further comprising the Vtep of: 
reconstructing said second pyramid after said uniform scaling step. 

25 4. The method of claim 2, wherein said determining a^et of optimal scale factors 
includes the steps of: 

determining a first voxel size in said first image; 
determining a second voxel size in said second image; and 
constructing a scaling matrix based on said first voxel siz^and said second 

30 voxel size. 

5. The method of claim 2, wherein said determining a set of optirr^l scale factors 
includes the steps of: 

segmenting a first anatomical structure in said first image; 
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segmenting a second anatomical structure in said second image; 
determining^ first largest span in said first anatomical structure; 

determining a second largest span in said second anatomical structure; and 

/"> 7 \ 

/ yC constructing a stealing matrix based on said first largest span and said second 



5 largest span. 

6. The method of claiik 2, wherein said determining a set of optimal translation 
and rotation parameters includes the steps of: 

segmenting a first anatomical structure in said first image to obtain a first map; 
10 segmenting a second anatomical structure in said second image to obtain a 

a second map; 

y2 identifying a first set of identifiable points or contours in said first map and a 

,j3 second set of identifiable points or contours in said second map; 

constructing a transformation matrix based on said first set of identifiable 
03 15 points or contours and said second set of\identifiable points or contours; and 
r* transforming said second image ^accordance with said transformation matrix 

to obtain a transformed second image. 



i 7. The method of claim 6, further comprising the step of: 

r* 20 subsampling said first image and said second image based on said first map 

and said second map. 

8. The method of claim 6, further comprising theVsteps of: 
measuring an error between said transformed second image and said first 

25 image; 

comparing said error to a threshold error; and 
repeating said identifying step if said error is greater Yhan said threshold error. 

9. The method of claim 6, further comprising the steps of: 
30 measuring an error between said transformed second ima^e and said first 

image; 

comparing said error to a previous error; and 
repeating said identifying step if said error is greater than said previous error. 
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10. The method ol\claim 1, wherein said locally transforming said first image and 

said second image includes the steps of: 

\ 

selecting a set of feature points in said first image; 

performing a feature tracking process on each of said set of feature points; and 
performing a transformation weighting process on each point in said first 
image that is not a feature point. 



11. A computer program product to use in conjunction with a computer system for 
registering multiple images, comprising: 

logic code for receiving a first image and a second image, said first image and 
said second image including three-dimensional data sets; 

logic code for globally transferring said second image via a motion tracking 
technique; 

logic code for locally transforming \aid first image via an iterative motion 
tracking technique; and 

logic code for outputting a transformed^first image and a transformed second 

image. 

12. The computer program product of claim 1 l\wherein said logic code for 
globally transforming includes: 

logic code for constructing a first pyramid for &aid first image; 

logic code for constructing a second pyramid fonsaid second image; 

logic code for determining a set of optimal scale factors based on said first 
pyramid and said second pyramid; 

logic code for performing a uniform scaling on said second image in 
accordance with said set of optimal scale factors to obtain a repealed second image; 

logic code for determining a set of optimal translation air^ rotation parameters; 

and 

logic code for performing a transformation on said rescaled^econd image 
based on said set of optimal translation and rotation. 

13. The computer program product of claim 12, further comprising l^gic code for 
reconstructing said second pyramid after said uniform scaling. 
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14. The computer program product of claim 12, wherein said logic code for 

determining a set of optical scale factors includes: 

logic code for determining a first voxel size in said first image; 
\ 

/ logic code for determining a second voxel size in said second image; and 

\ 

5 logic code for constructing a scaling matrix based on said first voxel size and 

said second voxel size. \^ 

15. The computer program product of claim 12, wherein said logic code for 
determining a set of optimal scale factors includes: 

10 logic code for segmenting a first anatomical structure in said first image; 

logic code for segmenting a seqpnd anatomical structure in said second image; 
logic code for determining a first largest span in said first anatomical structure; 
logic code for determining a seconfi largest span in said second anatomical 
structure; and 

15 logic code for constructing a scaling i^iatrix based on said first largest span and 

said second largest span. 



16. The computer program product of claim lv, wherein said logic code for 
determining a set of optimal translation and rotation parameters includes: 

20 logic code for segmenting a first anatomical structure in said first image to 

obtain a first map; 

logic code for segmenting a second anatomical Structure in said second image 
to obtain a second map; \ 

logic code for identifying a first set of identifiable points or contours in said 
25 first map and a second set of identifiable points or contours in said second map; 

logic code for constructing a transformation matrix baWd on said first set of 
identifiable points or contours and said second set of identifiable points or contours; 
and 

logic code for transforming said second image in accordance with said 
30 transformation matrix to obtain a transformed second image. 

17. The computer program product of claim 16, further comprising^ 
logic code for subsampling said first image and said second ima^ based on 

said first map and said second map. 
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18. The computer program product of claim 16, further comprising: 

logic code for measuring an error between said transformed second image and 

said first image; \ 

logic code for comparing said error to a threshold error; and 
logic code for repeating said identifying step if said error is greater than said 

threshold error. 



19. The computer program prodbct of claim 16, further comprising: 
logic code for measuring an error between said transformed second image and 

said first image; 

logic code for comparing said error to a previous error; and 
logic code for repeating said identifying step if said error is greater than said 

previous error. \^ 

20. The computer program product of claim 1\1, wherein said logic code for 
locally transforming said first image and said second image includes: 

logic code for selecting a set of feature points\in said first image; 
logic code for performing a feature tracking profess on each of said set of 
feature points; and 

logic code for performing a transformation weighting process on each point in 
said first image that is not a feature point. 
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